"Coulomb+linear" form of the static QCD potential in operator-product-expansion 
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The static QCD potential is analyzed in operator-product-expansion within potential-NRQCD 
framework when r <C Aq CD . We show that the leading short-distance contribution to the potential, 
defined as a perturbatively computable Wilson coefficient, can be expressed, up to 0(r 2 ), as a 
"Coulomb+linear" potential, ft coincides with the "Coulomb+linear" potential obtained previously 
from renormalon-dominance hypothesis. Non-perturbative contributions are 0(r 2 ) and subleading. 
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For decades, the static QCD potential Vqcd(V) has 
been widely studied for the purpose of elucidating the 
nature of the interaction between heavy quark and an- 
tiquark. It is defined from an expectation value of the 
Wilson loop as 
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where V is a rectangular loop of spatial extent r and time 
extent T. Generally, Vqcd(^) at short-distances can be 
computed accurately by perturbative QCD. On the other 
hand, the potential shape at long-distances should be de- 
termined by non-perturbative methods, such as lattice 
simulations or phenomenological potential-model analy- 
ses where phenomenological potentials are extracted from 
the experimental data for the heavy quarkonium spectra. 
Empirically it has been known that phenomenological po- 
tentials and lattice computations of Vqcd if) are both ap- 
proximated well by the sum of a Coulomb potential and 
a linear potential in the intermediate-distance range. 

Since the discovery of the cancellation of (D(Aqcd) 
renormalons in the total energy of a static quark- 
antiquark pair Etot{f) = Vqcd(?") + 2m po ie, convergence 
of the perturbative series for i^tot if) improved drastically 
and much more accurate perturbative predictions for the 
potential shape became available. It was understood that 
a large uncertainty originating from the O(Aqcd) renor- 
malon in Vqcd (t) can be absorbed into twice of the quark 
pole mass 2m po i c . Once this is achieved, a perturbative 
uncertainty of E to t(r) is estimated to be 0(Aq CD r 2 ) at 
r <C AqJ, d 0, based on the renormalon-dominance hy- 
pothesis. 

An operator-product-expansion (OPE) of Vqcd (r) was 
developed within an effective field theory (EFT) "po- 
tential non-relativistic QCD" (pNRQCD) Q. The idea 
of OPE is to factorize short-distance contributions into 
Wilson coefficients (perturbatively computable) and non- 
perturbative contributions into matrix elements of oper- 
ators, when the following hierarchy of scales exists: 



Aqcd <M/ « Vr. 
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Here, /i/ denotes the factorization scale. In this frame- 
work, residual renormalons, starting from 0(Aq CD r 2 ), 



are absorbed into the matrix elements of non-local 
operators (non-local gluon condensates). Then, in 
the multipole expansion at r <C AqJ, d , the leading 
non-perturbative contribution to the potential becomes 

0(A 3 qcd- 2 ) i 

Subsequently, several studies 0, Q showed that 

perturbative predictions for Vqcd(?") agree well with 
phonomenological potentials and lattice calculations of 
Vqcd(^), once the O(Aqcd) renormalon contained in 
Vqcd (r) is cancelled. Ref. 0] showed that a Borel resum- 
mation of the perturbative series gives a potential shape 
which agrees with lattice results, if the O(Aqcd) renor- 
malon is properly taken into account. In fact these agree- 
ments hold within the C(Aq CD r 2 ) uncertainty. These ob- 
servations support the validity of renormalon dominance 
and of OPE for FqcdO). 

Once the O(Aqcd) renormalon is cancelled, the per- 
turbative QCD potential becomes steeper than the 
Coulomb potential as r increases. This feature is un- 
derstood, within perturbative QCD, as an effect of the 
running of the strong coupling constant 0, Q . 

Moreover, using a scale- fixing prescription based on 
renormalon dominance hypothesis, it was shown ana- 
lytically that the perturbative QCD potential ap- 
proaches a "Coulomb+linear" form at large orders, up 
to an 0(Aq CD r 2 ) uncertainty. The "Coulomb+linear" 
potential can be computed systematically as more terms 
of perturbative series are included via renormalization- 
group (RG); up to the next-to-next-to- leading logarith- 
mic order (NNLL), it shows a convergence towards lattice 
results. 

In this paper, we analyze the QCD potential using 
OPE developed in j3[ and compare the leading Wilson 
coefficient (singlet potential) with the "Coulomb+linear" 
potential obtained in 0- 

The V-scheme coupling in momentum space ay (3) is 
defined as 
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where q = \q\; Cf is the second Casimir operator of 
the fundamental representation. In perturbative QCD, 
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av{q) is calculable in a series expansion of the strong 
coupling constant: 
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where, P n (£) denotes an n-th-degree polynomial of £ and 
o„ = _P fl (0). In this paper, unless the argument is spec- 
ified explicitly, as = as(/i) denotes the strong coupling 
constant renormalized at the renormalization scale /i, de- 
fined in the MS scheme. Here and hereafter, a v T (q) rep- 
resents a perturbative evaluation of ay(q) supplemented 
by RG evolution of as(q)- For instance, by a v T (q) up to 
NNLL, we mean that in © the sum is taken for n < 2 
and the three- loop running coupling is used for as(q)- 

The "Coulomb+linear" potential Vc+lM obtained in 
9], up to NNLL, is given by 
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where f3 n represents the (n+l)-loop coefficient of the beta 
function; e.g. in SU(3) Yang-Mills theory, (3q = 11, fo = 
102, • ■ •. The integral paths C\ and Ci are displayed 
in Fig. 2 of 0- The coefficient of the linear potential 
a can be expressed analytically in terms of the Lambda 
parameter in the MS-scheme A^jg. The "Coulomb" po- 
tential has a short-distance asymptotic behavior con- 



sistent with RG, Vc(r) 



-271-0 CM 



In 



r Ar- 



il 



In In 



, whereas its long-distance behavior 

is given by Vc(r) ~ —AttCf / (for); in the intermediate 
region both asymptotic forms are smoothly interpolated. 

Let us first present an intuitive argument. In fact, 
it already embraces an essential part of our discussion. 
We separate the integral (@J into the regions q > /if and 
q<lif. In the former region, ay (q) can be approximated 



well by a y T (q), 



hence we define 
2C F 



dq 



sin(qr) 
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Since /i/ ^> Aqcd, we expect that an accurate pertur- 
bative prediction for V\jy(r; /if) can be made. In the 
region q < /j,f, av(q) cannot be evaluated reliably in 
perturbation theory; rather it should be determined non- 
perturbatively. On the other hand, we may expand in r 
since fj.fr <C 1 by eq. J2J : 
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FIG. 1: Integral path C3 in the complex q-plane. q, denotes 
the IR singularity of as(q)- For 1-loop running, g» is a pole; 
beyond 1-loop running, g* is a branch point. In the latter 
case, branch cut is on the real axis starting from q» to —00. 



We can show that 

Vw(r;nf) - Vc+Ur) = const. + 0(fj, f r 2 ). (12) 

Eqs. jnj and l|12|l imply that the "Coulomb+linear" part 
of the QCD potential is determined by the short-distance 
contributions (q > (if), hence it is predictable in per- 
turbative QCD, while the non-perturbative contributions 
are of order jj^r 2 and subleading at r <C fij . [Through- 
out this paper, we are not concerned about the constant 
part of Vqcd(Oj keeping in mind that it can always be 
absorbed into 2m po i c in the total energy E to t (r) .} 
Eq. (|12fl can be shown as follows. According to JSJ), 
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where the integral path C3 is shown in Fig. ^ Since 
fifr 1, we may expand the Fourier factor as e %qr = 
1 + iqr — h (qr) 2 + ... in the above integral. Then one can 
show by suitable change of variables that the r -1 term 



2 ^lm/ dq^ 
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equals —A-kCf /(for) (at least) up to NNLL. Similarly, 
one can show for the r 1 term 
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(q). (10) Therefore, only remaining terms on the right-hand-side 
of (O are const. + 0(fj, 3 f r 2 ). 

Being intuitive, the above argument is subject to some 
flaws: (i) A factorization of scales is introduced only in 
the integral over q, whereas in a consistent OPE one 
should factorize scales in all quantum effects, namely in 
the computation of ay(q) as well, (ii) It is known that 
a n (n > 3) in a v T (q) includes IR divergences |lCj, so 
that Vtjv (r; fif) is not well-defined beyond NNLL. (iii) 
The perturbative series of V\jv(r; /jf) may still be an 
asymptotic series, hence one should clarify how to de- 
fine Vjj\r(r; (if). All these points (i)-(iii) are remedied in 



a consistent framework of OPE, in dimensional regular- 
ization and with appropriate renormalization procedure. 

An OPE of Vqcd (r) was developed in Q . In this and 
the next paragraph, we review the content of that pa- 
per relevant to our analysis. Within this framework, 
short-distance contributions are contained in the poten- 
tials, which are in fact the Wilson coefficients, while non- 
perturbative contributions are contained in the matrix 
elements that are organized in multipole expansion in f 
at r <C AqJ, d . The following relation was derived: 



V QCD (r) = V s (r)+SE vs (r) 

rri pOO 

SE vs = -ig 2 4^- I rfie-'^' 



(16) 



N c Jo 

+ 0(r 3 ). (17) 

Vs(r) denotes the singlet potential. i5£'us( r ) denotes 
the non-perturbative contribution to the QCD potential, 
which starts at 0(AQ CD r 2 ) in the multipole expansion. 
AV(r) = Vo{r) — Vs(r) denotes the difference between 
the octet and singlet potentials; see Q for details. In- 
tuitively Vs(r) corresponds to V\jv(r;nj) an d 8E\js(r) 
to Vm(r; fif). We adopt dimensional regularization in 
our analysis; we also refer to hard cutoff schemes when 
discussing conceptual aspects. 

In perturbative expansion in as, the QCD potential 
UqcdM coincides with the singlet potential Vs(r), i.e. 
SEtjs — 0. As already mentioned, perturbative expan- 
sion of Vqct>(t) in as includes IR divergences beyond 
0(a s ), hence Vs(r) also includes IR divergences in di- 
mensional regularization. SE\js(r) is expected to be non- 
zero beyond perturbation theory. In fact, if we do not 
expand AV(r) in a s in (|T7J) [la ] (but expand all other 
factors), 5E vs (r) becomes non-zero since AV(r) acts as 
an IR regulator. In this case, SE\js(r) contains UV di- 
vergences, given as poles in e, which exactly cancel the 
poles corresponding to the IR divergences in Vs (r) . Con- 
sequently, in the sum l|16[l. VqcdM becomes finite as 
e — > 0. These divergences in Vs(r) and 5Eus(r), respec- 
tively, can be regarded as artefacts of dimensional regu- 
larization, where the integral regions of virtual momenta 
extend from to oo. If we introduce a hard cutoff to each 
momentum integration, corresponding to the factoriza- 
tion scale /if, Vs(r) (q > (if) and 5E\js(r) (q < [if), re- 
spectively, would become finite and dependent on/if. In 
dimensional regularization, Vs(r) can be made finite by 
multiplicative renormalization, i.e. by adding a counter 
term (Z s - l)V s (r). 

With respect to the spirit of factorization in OPE, it 
is natural to subtract IR renormalons from Vg(r) in a 
similar manner. In 11], this was advocated and in prac- 
tice subtraction of the O(Aqcd) renormalon was car- 
ried out explicitly. The known IR renormalons of Vs(r) 
[=perturbative expansion of Vqcd(p)] are contained in 
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FIG. 2: Graph for SE vs {r) at 0(r 2 ) in the large-/3 ap- 
proximation. Bubble chain represents the renormalized gluon 



propagator in the large-/3o approximation, 
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In a hard cutoff scheme, it was shown 3] that the 
0(AQ CD r 2 ) IR renormalon of Vs(r) can be absorbed into 
5E\js(r). In dimensional regularization (D = 4 — 2e), 
one may compute, for instance, SEus(r) at 0{r 2 ) in the 
large-/?o approximation [13J, correponding to the graph 
in Fig. |21 It is given by 
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We note that the renormalon contribution (n! G n +\) and 
the UV divergences (multiple poles in e) are included 
in separate parts in this approximation. One can check 
explicitly that the 0(Aq CD r 2 ) UV renormalon, corre- 
sponding to the pole at u = 3/2 of G(u), cancels the 
0(Ag CD r 2 ) IR renormalon in Vs(r) |2| in the large-/?o ap- 
proximation. Therefore, in dimensional regularization, it 
is appropriate to subtract from Vs (r) the IR renormalons, 
e.g. in the form of l|18|) . 

We define a renormalized singlet potential (in dimen- 
sional regularization), in a scheme where the IR diver- 
gences and IR renormalons are subtracted, as 
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5av(q) is the counter term which subtracts the IR 
divergences, given as multiple poles in e, e.g. in the 
MS scheme. [l7j The dependence on (if is introduced 
through subtraction of the IR divergences and IR renor- 
malons. Then, we can apply our argument given through 
eqs. l|T5|l -(|15 fl to show that 

l/j R3 (r; M /) - V c+ L(r) = const. + 0{(i)r 2 ), (23) 

up to NNLL (since Sav(q) can be taken as zero up to this 
order). Moreover, beyond NNLL, this relation still holds 
after a simple replacement cty T (q) — > a^ T (g) + 5av(q) in 
the definition of Vc+h(r), © and ©• 

One may think that subtracting the integral (|18l) is not 
sufficient for subtracting all the IR renormalons. Our 
result l|23|) is unchanged, even if one subtracts the IR 
renormalon contributions using whatever other sophisti- 
cated method for estimating them. This is because the IR 
renormalons in Vs(r) take the form const. + 0(AQ CD r 2 ). 

The perturbative expansion of Vg R '(r; (if) may still be 
an asymptotic series. Since the IR renormalons have 
been subtracted and the factorization scale is set as 
(if ^> Aqcjq, we may expect that (r; (if) is Borel 
summable. 18] (At least, the Borel integral is convergent 
in the large-/3o approximation.) Then, we may define 
Vg (r) from the perturbative series either by Borel sum- 
mation or according to the prescription of |14I | ; both pre- 
scriptions lead to the same result when the series is Borel 
summable. 

Our result (|23[l shows that the renormalized singlet po- 
tential vj R) (r) can be expressed as a "Coulomb+linear" 
potential Vc+l(Oj u P to 0(/i^r 2 ), at short distances. We 
re-emphasize that there is no freedom to add a linear po- 
tential to Vc+lM in 123) 0- 

On the other hand, there is an arbitrariness in how 
to separate Vc+iS r ) into "Coulomb" and linear parts, 
as discussed in Stating more accurately, as yet we 
do not know any mathematically well-defined principle 
to separate Vc+lM into Coulombic and linear parts 
about r ~ 0, because of 1/lnr dependence in Vc+lM- 
Nonetheless, we consider the present separation Q © 
a natural one according to its construction, and also be- 
cause it is demonstrated that the perturbative prediction 
of Vqcd(^) up to O(ctg) is approximated well by this 



"Coulomb+linear" form for a fairly wide range of r and 

The author is grateful to A. Pineda and M. Tanabashi 
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